close all
clear all
clc

format longE

% Prima successione
% Inizializzazione di s a (vettore) zero di 100 elementi
s = zeros(1,100);
% caso base
s(1) = exp(-1);
% ciclo iterativo;
for i = 1:length(s)-1
    s(i+1) = 1-(i+1)*s(i);
end

% Seconda successione
% Inizializzazione di t a (vettore) zero di 100 elementi
t = zeros(1,100);
% caso base
t(length(t)) = 0;
% ciclo iterativo
for i = length(t)-1:-1:1
    t(i) = (1-t(i+1))/(i+1);
end


% Visualizzazione delle soluzioni
figure;
clf;
indvis = 1:30;
h = semilogy(indvis, abs(s(indvis)), 'k-');
set(h, 'LineWidth', 2);
hold on;
h = semilogy(indvis, abs(t(indvis)), 'm-');
set(h, 'LineWidth', 2);
hold off;
legend('s', 't');
title('Esempio di succ. da formulazione integrale');
xlabel('Iterazione');